Phase diagrams of Ising films with competing interactions 
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I. INTRODUCTION 

In recent years magnetism in thin films of a few atomic 
layers has attracted much interest, both theoretically and 
experimentallyEI. However, studies on the influence of 
the layer thickness, L, on spatialbj jnodulated magnetic 
structures seem to be very scarceo'El, albeit the possible 
lack of compatibility between the film geometry and the 
modulations in the bulk as well as the effect of the sur- 
faces on the ordering phenomena may lead to interesting 
features. 

In this article, we shall deal with this topic by analysing 
phase diagrams of the axial next -nearest-neighbour Ising 
(ANNNI) modeloQ on a simple cubic lattice. Due to its 
competing interactions, the model displays, in the limit 
of infinitely many layers, a phase diagram with a plen- 
itude of commensurate phases, including those spring- 
ing from the multiphase point at zero temperature and 
those emerging from structure combination branching 
processes at finite temperatures, as well as incommen- 
surate phases and a Lifshitz point. Many of these as- 
pects have been observed experimentally, in particular 
in magaete, alloys, polymers, and iepnelectrics, see the 
review^ H and related recent workCrllil. 

The intriguing structural complexities are severely af- 
fected when the lattice consists of rather few layers per- 
pendicular to the axis of competition. Of course, in the 
case of periodic boundary conditions for the top and bot- 
tom layerjS, jtbe bulk phases which fit to the film thickness 
still exista^BO. In the case of free boundary conditions 
for the surface layers, describing more realistically exper- 
imental situations, novel surface-induced features may 
evolve. The aim of this paper is to identify similarities 
and typical differences between the phase diagrams in 
the limit of an infinite lattice and for thin films with free 
boundaries. 

To study the influence of the surfaces in more detail, 
we also varied the intralayer couplings in the surfaces 
compared to those in the bulk. In semi-infinite systems, 
one then encounters the well-known surface critical phe- 
nomena with eithp; pcdinary or extraordinary and surface 
phase transitional3'Eil, depending on the strength of the 
surface and bulk couplings. Note that critical properties 
of the surface magnetization at the Lifshitz point in the 



semi-infinite ANNNI model have beep_analysed recently, 
refining results o£-rnean-field theoryllj by doing Monte 
Carlo simulational3. 

In the following, we shall mostly deal with thin films 
of up to ten layers. For each film thickness, dis- 
tinct phase diagrams in the (temperature-competition 
strength)-plane are determined, using mean-field theory, 
low-temperature expansions, and Monte Carlo simula- 
tions. Mean-field theory is found to usually provide reli- 
able guidance to the correct phase diagrams; qualitative 
shortcomings are observed especially in the case of van- 
ishing surface couplings. Brief accounts of some of our 
findings for le^Pfl surface and bulk couplings have been 
given beforetZmS. 

The article is organized as follows: first the model and 
the methods are presented, then the resulting phase dia- 
grams are discussed. The paper is concluded by a sum- 
mary. 



II. MODEL AND METHODS 

We consider the ANNNI model on a simple cubic lat- 
tice (setting the lattice constant equal to one) for films 
of L layers, L > 2, with free boundary conditions for 
the surface layers. Each spin Sj, at site j, can take 



only two values, Sj — 1 (spin 'up') or Sj 



-1 (spin 



'down'). The interactions are supposed to be ferromag- 
netic between neighbouring spins in each layer, Jq > 0, 
as well as between neighbouring spins in adjacent layers, 
Ji > 0, and to be antiferromagnetic, J2 < 0, between 
axial next-nearest neighbour spins, distinguishing one of 
the three cubic axes, say, the ^-axis. The strength of the 
competition between the interactions along the z~axis is 
K — — J2/J1. The intralayer couplings Jq may be differ- 
ent in the (top and bottom) surface layers, denoted by 
Js, compared to those in the bulk layers, Jb, with the ra- 
tio r — Js/ Jb- For simplicity, we set Jb = Ji- Changing 
the interlayer coupling Ji relative to the bulk intralayer 
interaction, Jb, modifies onljz—quantitatively the phase 
diagram in the limit L — > oot£l"EJ. 

The ground state properties of the ANNNI films are 
readily obtained. We first consider Jg > 0. The spins in 
all layers are aligned ferromagnetically at k < 1/2. For 
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K > 1/2, and even film thickness L, the ground state 
consists otflairs of layers with, say 'up' spins, denoted as 
a '2-band'lH3'E2l, followed by a 2-band of 'down' spins, be- 
ing obviously two-fold degenerate by interchanging the + 
and — spins. If L is odd, then for 1 > k > 1/2, the ground 
state comprises one 3-band and (L — 3)/2 2-bands, being 
(L — l)-fold degenerate. For k > 1, structures with one 
1-band at the surface and 2-bands are stable, with a 4- 
fold degeneracy. Both for even and odd number of layers 
L, the degeneracy Dl of the multiphase point at k = 1/2 
is given by a Fibonacci sequence Dl = Dl-i + 
where D2 = Z?3 = 2, corresponding to all possible combi- 
nations of A;— bands with fc > 1, as discussed beforecj'EJ. 
For L odd, at k = 1, L + 3 structures have the same 
ground state energy. For — Q, the situation at k > 1 
is of special interest, L being odd. At k = 1, structures 
with 2-bands followed by an arbitrarily ordered surface 
layer are degenerate, due to the compensation of the in- 
terlayer interactions between the surface spins and those 
in the adjacent and next-nearest layers. In the case of 
L = 3 and k > 1, the orientations of the spins in one of 
the surface layers are random, and completey antiparal- 
lel to that random pattern in the other surface; in the 
center layer, all spins have the same sign. 

At non-zero temperatures, the stable phases have been 
determined using mean-field theory, low-temperature 
expansions, especially about the special points (T = 0, 
K = 1/2) and (T = 0, k = 1), as well as Monte Carlo 
techniques. 

In mean-field theory, the magnetization per layer, m^. 



i = 1,2,. ..L, in 
dard equations!!. 



. equilibrium follows from the stan- 



rui = tanh((4Jo(i)mi + Ji(mj_i + mi+i) 

+J2{mi-2 + ni.,+2))/kBT) (1) 

where Jo(i) = Js for z = 1, L and Jo(?) — Jb for the other 
layers. Free boundary conditions are implemented by set- 
ting rrii = for i = —1, 0, L -f- 1, L + 2. To obtain the 
thermally stable magnetization pattern, we solved Eq. 
(1) iteratively starting from all possible combinations of 
fully ordered, = ±1, or completely disordered layers, 
rrii = 0, i.e. from, in principle, 3^ distinct configurations 
(the number may be reduced using symmetry consider- 
ations). We then determined among the solutions the 
one with the lowest free energy. Phase boundaries are 
identified by singularities in the free energy and specific 
heat. Obviously, a fine scan of the (temperature ksT / Ji, 
competition strength K)-plane is rather computer-time 
consuming, and the full analysis of the mean-field the- 
ory was done for films of up to L = 10 layers. The ratio 
r = Ja/ Jb varied from to 1.5, which would cover, in the 
limit L — > 00, both ordinary and surface transitions. 

The transition to the paramagnetic phase may be stud- 
ied for films of larger thickness, L, by analysing the lin- 
earized form of the mean-field theory. From Eq. (1) 
one obtains a matrix of rank L, with the eigenvalues de- 
termining the phase transition temperature, Tc, and the 



eigenvectors describing the critical magnetization pat- 
tern. Films of thickness up to L = 50 were considered, 
especially for competition ratios close to that of the Lif- 
shitz point in the infinite system, i. e. k = 0.25. 

The phase diagram close to the special ground states 
at K = 1/2 and 1 may be iavestigated by using exact 
low-temperature expansionsE3E3. The stable phases are 
identified by calculating the free energy resulting from 
spin excitations for all structures being degenerate at the 
two special points. Indeed, to establish the stable phases 
springing from the special points for thin films (we stud- 
ied films with up to ten layers), it usually suffices to do 
expansions up to first order, involving merely a single 
spin flip. 

Complementary to the low-temperature analysis, 
Monte Carlo simulations may be applied to study the 
phase diagrams at higher temperatures. We used both 
the standard single-flip Metropolis algorithmn3 and a 
cluster-flip algorithmo (attention may be drawn to an- 
other cluster-flip algorithm to simulate spatially mod- 
ulated structures in Ising modelstj). The main aim of 
our simulations has been to check results of mean-field 
theory, restricting ourselves to selected cases. Of course, 
the layers are now finite, consisting of, say, spins. 
Thereby, finite-size effects have to be taken into account 
(here, attention may be drawn to early Monte Carlo work 
on the nearest neighbour Ising model on thin filmsEil). 
Typically M was varied from 10 to 100. To identify the 
structures and boundary lines of the various phases, we 
computed the energy, the specific heat, the layer magne- 
tizations and corresponding histograms as well as corre- 
lation functions between spins in different layers. Each 
run was performed with about 10^ Monte Carlo steps 
per spin, when using the single-spin-flip algorithm. In 
the cluster-algorithm, in each run usually about 2 • lO'^ 
clusters were generated. 



III. RESULTS 

For each film thickness L, distinct phase diagrams are 
obtained. Before presenting them, we shall outline some 
general features, as inferred from calculations for films 
with up to ten layers. 

One may distinguish between phases corresponding to 
ground states and those which are stabilized only at 
higher temperatures. Especially, the ferromagnetic phase 
as well as phases consisting of 2-bands augmented, for 
L odd, by one 3-band or one 1-band, evolve from the 
ground states. A 'A:-band', at T > 0, means that spins in 
k adjacent layers are oriented predominantly in one direc- 
tion, preceded and followed by layers where the spins are 
oriented predominantly in the opposite direction. The 
3-band tends to be at or near to the center of the film at 
low temperatures, when it gains entropy from easy spin 
flips in the two border layers of the 3-band. The band, 
however, sticks to one of the surfaces of the fllm when J, 
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is rather weak so that entropy is gained from energeti- 
cally easy flips of surface spins, as may be readily seen 
from low-temperature expansions. 

Phases corresponding to ground states with 2-bands 
and one 3-band at or near the center (or at the surface) 
of the film may be stabilized only at higher temperatures, 
due to entropic or symmetrization reasons, while at low 
temperatures the 3-band sticks to the surface (or center) 
of the film. 

The phases with one 3-band and one 1-band, L odd, 
are separated, for Jg > 0, by the transition line arising 
from (r = 0, K = 1); the line is at low temperatures of 
first order. If the surface couplings Jg are not very strong, 
the line terminates in a critical point, above which the 
two phases become indistinguishable. That point moves 
down to zero temperature as Js vanishes. 

Other stable low-temperature phases may spring from 
the multiphase point (T = 0, k = 1/2) for L > 5. In 
fact, in the limit L — > cx3, a sequence of infinitely 
many periodic (32*) phases, i =0,1,2,. ..,cx3, arises from 
that point (here and in the following, we are using or 
adapting the standard notationcHlc3, setting the sequence 
of bands in a phase in '() '-brackets). For thin films, 
the low-temperature expansions indicate the following 
systematics: If L is a multiple of 3, the (3'^/'^) struc- 
tures become stable next to the ferromagnetic phase; for 
L = 3n -I- 1(= 3n -I- 2), n > 1, they are replaced by 
the (3"^^4) ((3"^^5)) structures. Further phases arising 
from the multiphase point show up for L > 7, having, 
for L — 8, two 3~bands and two 2-bands, and two 3- 
bands and three 2-bands for L = 10. The 3-bands stick 
to the surfaces of the film, when is sufficiently weak, 
as before. An increasing number of similar phases com- 
prising 2-bands and 3~bands are expected to occur for 
thicker films. In comparison to the (32*) phases, arising 
from the multiphase point in the limit L — > oo, the se- 
quences of 2-bands and 3-bands are slighty rearranged 
and modified, refiecting the influence of the surfaces. 

The transitions between the phases emerging from the 
multiphase point are, at low temperatures, of first order. 
Note that, for films, some, but not all these phases extend 
up to the transition line to the paramagnetic phase, Tc- 

Novel phases may become stable close to Tc- As may 
be obtained readily from the linearized mean-field the- 
ory, confirmed by Monte Carlo simulations, the ordered 
phases directly below this transition line are alternat- 
ingly, as k increases, symmetric and antisymmetric with 
respect to the center plane (being obviously a real layer 
for L odd, and a fictitious plane for L even), changing 
thereby the parity. Each change is associated with an 
abrupt decrease in the average wavelength, as follows 
from a Fourier analysis of the magnetization pattern, m^. 
For odd L, the magnetization in the center layer vanishes 
in the antisymmetric cases, denoted as 'Ig-band'. There, 
the interlayer couplings with the neighbouring layers tend 
to compensate. Such 'partially disordered phases' are, of 
course not stable below the ordering temperature of the 
two-dimensional Ising plane. They do not exist in the 



limit L — > ooEaEil.- Of course, as usual, the location of 
Tc is largely overestimated in mean-field theory, as seen 
from comparison with simulational data. 

An interesting feature of the transition line to the 
paramagnetic phase is the point where the ferromagnetic 
structure becomes unstable against spatially modulated 
structures, i.e. the Lifshitz point in the limit L — > oo. In 
mean-field theory, for rather thick films of up to 50 lay- 
ers, the surface magnetization is found to show intriguing 
crossover effects with an effective critical exponent close 
to 1 in the vicinity of that point, in agreement with the 
mean-field analysis on the semi-infinite ANNNI modcllla. 
In a recent Monte Carlo simulation, the exponent has 
been shown to be significantly affectedphy critical fluctu- 
ations, acquiring a value of about 0.62ll3. 

Further phases may show up at intermediate temper- 
atures, in particular, the (54) and (64) structues for 
L = 9 and 10. They are reminiscent of the plenitude 
of commensurate phases on the ferromagnetie side of the 
ANNNI model in the thermodynamic limitcZl, and more 
of them are expected to occur in thicker films. Likewise, 
wider films are needed to see indications of the systematic 
branching processes involving structures with 2-bands 
and 3-bands, as occurs for L — > ooEZl. 

Note that the magnetization in a given phase may 
change sign in some layers when varying temperature 
and/or strength of couplings, keeping the symmetry of 
the magnetization pattern about the center plane (sym- 
metric, antisymmetric or asymmetric). For instance, for 
an (anti-)symmetric phase such a reversal of the sign 
may happen in the magnetization of the surface layers 
belonging originally to fc— bands with fc > 2, when k is 
increased, especially for weak surface couplings. In any 
event, we characterise a phase by the sequence of bands 
which occurs at lowest temperatures in case of equal sur- 
face and bulk couplings, following the unambiguous no- 
tation for the infinite latticec^. 

We shall now discuss results for each film thickness, 
from L = 3 to L = 10, illustrating the general features 
in selected cases in more detail. 

In films with three layers, L — 3, the ferromagnetic, 
the (21) and the (HqI) phases show up. As found in 
simulations, for equal surface Jg and bulk Jj, couplings, 
the ferromagnetic phase is stable at k < 1, the other two 
phases occur at k > 1, with the transition between the 
(21) and the partially disordered phases being of second 
order, see Fig. 1, and belonging to the universality class 
of the two-dimensional Ising model. At k ~ 1, the fer- 
romagnetic and (21) phases coexist, with the transition 
being of first order. By weakening the surface interac- 
tions, Js, the ferromagnetic phase extends to larger val- 
ues of K, due to the entropic effect in the surface layers. 
In the limit Jg — 0, the (21) and (HqI) phases are com- 
pletely suppressed (this aspect is not correctly described 
in mean-field theory). On the other hand, by enhancing 
the surface couphngs, the (21) and the (HqI) phases are 
favoured and stabilized already at k < 1, with the tran- 
sitions to the ferromagnetic phase being always of first 
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order. 
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FIG. 1. Phase diagram of the ANNNI model for L = 3 
and Js — Jb using Monte Carlo simulations. The curves are 
guides to the eye. Phase boundaries of first (second) order 
are denoted by dashed (solid) lines. 

For L = A, the phase diagrams are very simple, with 
two ordered phases, the symmetric, ferromagnetic and 
the antisymmetric (2^) phases being separated by a first 
order transition fine arising from the special point (T — 
0,K = 1/2). By weakening Js, the range of stability of the 
symmetric phase is again enlargened. The (2^) structure 
is obviously favoured when the surface layers are more 
ordered, and thence its range of stability expands when 
Js gets stronger. 

Somewhat richer phase diagrams are encountered when 
L = 5. For Js = Jb, one finds, increasing the competition 
ratio K, at low temperatures the ferromagnetic, the (32), 
and the (2^1) phases, with the latter two becoming unsta- 
ble against the (21o2) and (131) phases at higher temper- 
atures. Phase diagrams, obtained from mean-field theory 
and Monte Carlo simulations have been depicted before 
for the case Js — JblUl. The transition line between the 
two asymmetric (with respect to the center plane) (32) 
and (2^1) phases ends in a critical point above which the 
phases become essentially indistinguishable. By weaken- 
ing Js, the critical point moves to lower temperatures. 
Eventually, in the limit of zero surface couplings, there is 
only one asymmetric phase, with one of the surface mag- 
netizations going to zero on approach to (T = 0, k = 1). 
The phase diagram, as determined from simulations, is 
shown in Fig. 2a. Note that the (131) structure, be- 
ing present close to Tc in the mean-field calculations, is 
squeezed out (as we checked for k as large as 15). In 
turn, the range of stability of the (21o2) phase is under- 
estimated in mean-field theory, similar to the situation 
for equal surface and bulk couplingalZl. The magnetic 
disorder in the center layer sets in roughly at the tran- 
sition temperature of the two-dimensional Ising model. 



because then the interlayer couplings Ji and J2 to the 
axial spins in the two adjacent and two next-nearest lay- 
ers are largely cancelled. When increasing the surface 
couplings Js, the (131) phase is stabilized, as shown in 
Fig 2b, in qualitative agreement with mean-field theory. 
Actually, the magnetization then tends to be close to zero 
in the second and fourth layers, being rather large with 
the same sign in the surface layers and with the opposite 
sign in the center layer. 
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FIG. 2. Simulated phase diagram of the ANNNI model for 
films with five layers at Js/ Jb= (a) and (b) 1.5. 

For L = 6, only three ordered structures exist, the 
ferromagnetic, the (3^), and the (2'^) phases. The tran- 
sitions between them are of first order. By decreasing 
the surface interactions, the range of stability of the (3^) 
phase gets wider, with 1-bands at the surfaces for large 
competition ratio k. This range and that of the ferromag- 
netic phase shrink when enhancing the surface couplings, 
and the (2^^) structure is favoured. Predictions of mean- 
field theory were qualitatively confirmed in simulations. 
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Adding another layer to the film, L = 7, leads to quite 
different phase diagrams, with possibly seven distinct or- 
dered phases, namely the ferromagnetic, (43), (31o3), 
(232), (322), ^231), and (121o21) phases, as found in 
mean-field theory. For Jg — Jb, the (232) phase extends 
down to zero temperature at 1/2 < k < 1, while the 
(32^) structure is stabilized only at intermediate temper- 
atures, where the entropic effect of-the small surface mag- 
netization plays the crucial rolellll. Examples of phase 
diagrams for different surface and bulk couplings are de- 
picted in Figs. 3a and 3b, for r = Js/Jb = 0.25 and 
1.5, as obtained from mean-field theory. Compared to 
the case of Js = Jb, the enhancement or weakening of 
the surface couplings may lead to similar trends as for 
i = 5, but also to novel features. As before, when 
Js— 0, there is no transition line arising from the spe- 
cial point (T = 0, K = 1), i.e. the two asymmetric (32^) 
and (2'^1) structures transform gradually into each other 
at non-zero temperatures, with a vanishing magnetiza- 
tion in one of the surface layers on approach to the spe- 
cial point. Indeed, even the transition line to the neigh- 
bouring asymmetric (43) phase terminates in a critical 
point, above which all these asymmetric structures be- 
long to the same phase. Perhaps somewhat unexpect- 
edly, at 1/2 < K < 1, the (2^3) structure forms the 
low-temperature phase when r < (3 -I- k)/4, while for 
stronger surface couplings the symmetric (232) phase is 
stable down to zero temperature. In the former case, the 
symmetric phase is stable next to Tc, see Fig. 3a. A sim- 
ilar entropy-driven shift of the 3-band from the interior 
of the film to one of its surfaces occurs at wider films, L 
odd, as well, for sufficiently weak Jg. The effect follows 
from the low-temperature analysis, and it is described 
correctly by mean-field theory. 

In the case of L = 8, six ordered phases are observed; 
the phase diagram as obtained from mean-field theocy 
for equal couplings Js ~ Jb has been shown beforeEII. 
Apart from the trivial ferromagnetic and (2**) phases, 
there is the (3^2) phase, which may, at sufficiently strong 
surface couplings Js, arise from the multiphase point 
(T = 0, K = 1/2), with the symmetric (323) struc- 
ture, formed at higher temperatures, remaining stable 
up to the transition line to the paramagnetic phase, Tc- 
The latter phase is stable down to the multiphase point 
for sufficiently weak surface couplings, where entropy 
is gained, again, from 3-bands at the surfaces, as fol- 
lows from low-temperature considerations, in agreement 
with mean-field calculations. Actually, for small Js, the 
(323) phase persists even for quite large values of the 
competition ratio k. Furthermore, the (53) phase may 
spring from the multiphase point. The antisymmetric 
(4^) phase is stable close to Tc, extending down to lower 
temperatures with decreasing surface couplings. 
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FIG. 3. Phase diagram of the ANNNI model for L = 7 
at Js/ Jb= (a) 0.25 and (b) 1.5, as obtained from mean-field 
theory. The critical point at the end of a line of transitions 
of first order is denoted by an open circle. 

In the case of nine layers, the rather complex phase 
diagrams consist of up to nine ordered phases. Close 
to Tc, increasing n, the alternatingly symmetric and an- 
tisymmetric phases are the ferromagnetic, (41o4), (3'^), 
(2^102^), and (12321) phases. At low temperatures, one 
may, in addition, encounter, the (2^1) phase as well as, 
for fairly weak surface couplings, r < (3-|-k)/4, the (32^) 
phase, and, for stronger couplings, the (2^32) phase at 
1/2 < ft < 1, similar to the situation in the case L = 7. 
Again, this feature follows from the low-temperature 
analysis, and it is described correctly by mean-field the- 
ory, see. Fig. 4 (the case Js = Jb has been depicted 
beforeE£l; there, the (32^) phase is squeezed out). Thence, 
by lowering the surface couplings, the 3-band eventually 
sticks at the surface. In fact, the same phenomenon holds 
for wider odd films as well. The transition line, sepa- 
rating the structures with one 1-band from that with 
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one 3-band and arising from (T = 0, k = 1), terminates 
again in a critical point for moderate and weak surface 
couplings. The point moves to zero temperature as J, 
vanishes. Perhaps most interestingly, mean-field theory 
suggests that the (54) phase is stable at intermediate 
temperatures, as shown in Fig. 4. Indeed, its existence 
has been confirmed in Monte Carlo simulations for equal 
surface, Js, and bulk, J^, couplings. The phase shows up 
below the (41o4} phase which, in turn, becomes stable at 
about the transition temperature of the two-dimensional 
Ising model, like the corresponding partially disordered 
phases for films with an odd number of layers. 
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FIG. 4. Phase diagram of the ANNNI model for L = 9 at 
Js/Jb~ 0.95 using mean-field theory. 

Finally, for L — 10, seven distinct ordered phases have 
been identified and located applying mean-field theory 
for Js = Jf,t3. Close to the transition to the paramag- 
netic phase, Tc, the ferromagnetic, (5^), (343), (23^2), 
and (2^) phases occur. The (43^) phase springs from 
the multiphase point (T = 0, k = 1/2). At intermediate 
temperatures, the (64) phase has been detected in a tiny 
region of the phase diagram. The phase persists, when 
varying the surface couplings. For sufficiently small val- 
ues of Js, the (32^3) phase is stabilized, at the cost of 
the (23^2) phase, near the multiphase point, reflecting 
again the gain of entropy due to easy excitations of sur- 
face spins belonging to 3-bands for structures comprising 
3-bands and 2~bands. 



IV. SUMMARY 

The ANNNI model on thin films has been studied, 
using mean-field theory, low-temperature analyses, and 
Monte Carlo techniques. For each film thickness, vary- 
ing the number of layers from L = 3 to L = 10, dis- 
tinct phase diagrams have been determined, monitoring 
surface-induced features by also changing the strength of 



the couplings in the surface layers. 

In the limit of infinite lattices, L — > oo, the com- 
peting interactions of the ANNNI model lead to a phase 
diagram with a rich variety of spatially modulated mag- 
netic structures. Signatures of some of these features, like 
the sequence of commensurate phases springing from the 
multiphase point, are already present in thin films, with 
the films displaying generic and distinct features as well. 
In particular: (i) The ordered phases occuring directly 
below the transition line to the paramagnetic phase are 
alternatingly, as the ratio n — — J2/J1 between the com- 
peting interactions of the model increases, symmetric and 
antisymmetric about the center of the film, leading to 
partially disordered phases with a paramagnetic center 
plane for films with an odd number of layers, (ii) The 
phases springing from the multiphase point may include 
structures consisting not only of 2-bands and 3-bands, 
but also with a 4-band or a 5-band. (iii) For L odd, 
there is an additional transition line near k ~ I, sepa- 
rating phases with 2-bands augmented by one 3-band or 
one 1-band. 

By modifying the surface couplings, additional inter- 
esting phenomena are induced. Indeed, the range of sta- 
bility of the various phases as well as the type of phase 
transition may be affected by varying the couplings. En- 
tire transition lines may, in fact, disappear. Perhaps most 
noticeable, 3-bands may be forced to stick to the surface 
or move to the interior of the film, thereby possibly form- 
ing new phases with a different symmetry. 

These findings may encourage experimental work on 
thin films of magnets or alloys showing complicated spa- 
tial orderings in the bulk. On the theoretical side, several 
extensions of the present study seem to be promising, in- 
cluding thejCffects of a larger number of layers, modified 
interactionsE3E3, and magnetic bulk as well as surface 
fields, in order to possibly detect additional systematics 
in the intriguing phase diagrams of the ANNNI model 
and variants. 
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